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IS.  ABSTRACT 

Project  Secede  is  an  ARPA  program  aimed  at  solving  certain  defense  related  pro¬ 
blems  through  the  study  of  chemical  releases  in  the  ionosphere.  The  chemical  of  pri¬ 
mary  interest  has  been  barium.  The  aspects  of  the  release  which  are  of  most  interest 
include  the  growth  of  the  cloud  and  the  evolution  of  its  striated  structure.  Several 
test  series  of  releases  have  been  conducted  the  latest  of  which  was  Secede  II. 

The  present  report  concerns  data  obtained  during  Secede  II  by  an  optical  inter¬ 
ferometer  built  and  operated  by  M.  A.  biondi,  R.  Hake,  and  D.  Sipler  of  the  University 
of  Pittsburgh  with  the  support  of  R.  Lambert  and  the  General  Electric  Space  Sciences 
Laboratory.  The  report  consists  of  a  number  of  data  volumes  and  this  one  summary 
textual  volume.  Only  the  summary  volume  will  receive  general  circulation;  the  data 
volumes  will  be  distributed  to  only  a  few  users  and  additional  copies  may  be  requeste< 
from  the  contracting  officer.  Data  volumes  exist  for  events  Spruce,  Olive,  Redwood, 
and  Plum.  > 

Interferometer  data  on  events  Spruce,  Olive,  Redwood  and  Plum  has  been  reduced 
and  photographs  of  representative  10  records  on  Spruce  have  been  made  10  records  of 
other  events  are  absent  presumably  because  of  non-optimum  adjustment  of  the  system 
during  those  events.  The  interferometer  data  is  presented  in  detail  in  the  Event 
Data  Books  published  separately  while  the  10  records  on  Spruce  are  reproduced  in  this 
volume. 
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I.  BACKGROLND  AND  PERSPECTIVES 


Project  Secede  is  an  ARPA  program  aimed  at  solving  certain  defense  related 
problems  through  the  study  of  chemical  releases  in  the  ionosphere.  The  chemical  of 
primary  interest  has  been  barium.  The  aspects  of  the  release  which  are  of  most  interest 
include  the  growth  of  the  cloud  and  the  evolution  of  its  striated  structure.  Several  test 
series  of  releases  have  been  conducted  the  latest  of  which  was  Secede  U. 

The  present  report  concerns  data  obtained  during  Secede  n  by  an  optical 
interferometer  built  and  operated  by  M.  A.  Biondi,  R.  Hake,  and  D.  Sipler  of  the 
University  of  Pittsburgh  with  the  support  of  R.  Lambert  and  the  General  Electric  Space 
Sciences  Laboratory.  The  report  consists  of  a  number  of  data  volumes  and  this  one 

summary  textual  volume.  Only  the  summary  volume  will  receive  general  circulation; 

» 

the  data  volumes  will  be  distributed  to  only  a  few  users  and  additional  copies  may  be 
requested  from  the  contracting  officer.  Data  volumes  exist  for  events  Spruce,  Olive, 
Redwood,  and  Plum. 

Briefly,  the  series  was  held  in  the  vicinity  of  Eglin  Air  Force  Base, 

Florida,  during  January  and  early  February  1971.  The  test  matrix  is  shown  in  Table  1. 


TABLE  1 

SECEDE  H  TEST  MATRIX  (U) 


PAYLOAD 


EVENT 

DATE 

RELEASE  TIME 
(CST) 

ALTITUDE 

(km) 

WEIGHT 

(kg) 

NUTMEG 

16  Jan  71 

1734:40 

150 

48 

OLIVE 

29  Jan  71 

1753:57 

185 

352 

PLUM 

20  Jan  71 

1747:06 

185 

48 

QUINCE 

2  Feb  71 

1034:06 

185 

48 

REDWOOD 

2r  71 

1752:08 

250 

48 

SPRUCE 

1  Feb  71 

1752:04 

185 

48 
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The  instrument  complex  has  been  described  in  detail  in  Reference  1  and 
2,  but  the  salient  characteristics  are  summarize'*  as  follows.  Instrumentation  con¬ 
sisted  of  three  major  devices :  a  two  cl  rmel  Fabry-Perot  Interferometer  (FP);  a  filter  photo¬ 
meter;  and  an  image  orthicon  system.  The  FP  was  used  to  obtain  spectral  line  profiles 
by  scanning  wavelength  (by  changing  driving  gas  pressure)  while  viewing  in  a  fixed 
direction  and  to  obtain  spatial  intensity  profiles  by  scanning  the  cloud  while  maintaining 
a  constant  wavelength.  The  filter  photometer  was  used  to  obtain  light  intensity  over  a 
relatively  wide  field  of  view  integrated  over  an  entire  spectral  line.  Normally  the  photo¬ 
meter  viewed  the  4934A  feature  but  had  the  capability  of  viewing  the  5535,  6111,  and 
6142  lines  as  well.  The  10  system  is  described  in  Reference  2. 

IL  DATA  REDUCTION 

Interferometer  data  on  events  Spruce,  Olive,  Redwood  and  Plum  has  been 
reduced  and  photographs  of  representative  IO  records  on  Spruce  have  been  made.  IO 
records  of  other  events  are  absent  presumably  because  of  non-optimum  adjustment  of 
the  system  during  those  events.  The  interferometer  data  is  presented  in  detail  in  the 
Event  Data  Books  published  separately  while  the  IO  records  on  Spruce  are  reproduced 
in  the  next  section  of  this  volume. 

The  first  element  of  the  reduced  interferometer  data  is  a  plot  of  azimuth 
and  elevation  tor  each  scan.  Instead  of  merely  presenting  a  rectilinear  presentation 
it  was  decided  that  the  pointing  information  could  better  be  arranged  as  a  scan  line 
which  could  be  laid  over  a  photograph  of  the  cloud  so  to  facilitate  identification  of 
structural  features  and  comparison  of  digital  records  and  densitometric  results.  A 
photograph  is  a  representation  of  the  projection  of  a  3-d  scene  on  tl  plane  (assuming 
no  distortion)  normal  to  the  optical  axis  of  the  system.  To  compa  e  the  interferometer 
scan  line  with  such  a  photograph  it  is  necessary:  a)  to  have  had  the  camera  and  inter¬ 
ferometer  sufficiently  close  to  each  other  that  they  viewed  the  cloud  from  nearly  the 
same  aspect  and  b)  to  present  the  scan  line  as  a  projection  upon  the  same  plane  that 
is  represented  by  the  photograph.  To  accomplish  this  the  scan  line  is  presented  in 
the  gnomonic  projection. 
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The  gnomonic  projection  is  the  projection  of  a  sphere  upon  a  plane  tangent 
to  the  sphere  at  a  given  point,  as  seen  from  the  center  of  the  sphere.  An  important 
property  of  this  projection  is  the  fact  that  all  great  circles  are  represented  by  straight 
lines. 

In  Figure  1  let  the  circle  represent  the  meridional  plane  of  the  sphere 
(i.  e.  constant  azimuth)  on  which  the  projection  plane  is  tangent.  Let  the  radius  of  the 
sphere  be  R.  The  point  A,  located  an  angular  distance  a  from  the  radius  vector  to 
the  projection  plane.  A*  is  located  at  a  distance  /  =  R  tan  a  from  the  point  of 
tangency  P.  All  points  located  at  an  angle  of  a  from  the  radius  OP  project  onto  a 
circle  of  radius  R  tan  a  on  the  projection  plane,  as  shown  in  Figure  2,  In  that  figure 
the  Y  axis  is  taken  to  be  the  projection  of  the  meridian  through  the  point  of  tangency 
while  the  X  axis  is  the  projection  of  the  orthogonal  great  circle  going  through  the  same 
point.  If  we  let 

4  *  elevation  of  the  mount  <with  the  equator  identical  with  horizon) 

A  *  azimuth  measured  from  the  meridian  of  tangency. 
and  the  elevation  of  the  pcint  of  tangency 

Then  a  point  located  at  {$,  X )  on  the  sphere  will  be  mapped  into  a  point  located  at 
(x,  y)  on  the  plane  of  Figure  2  where 

X  ^  Sin  X  Cos  $ _ _ 

R  *  Sin  4>0  Sin  4  +  Cos  Cos  A  Cos  $ 

and  Y_  m  Cos  Sin  $  -  Sin  Cos  A  Cos  $ 

R  Sin  $  Sin  $>0  +  Cos  $  Cos  $0  Cos  X 

Using  these  equations,  points  on  a  sphere  can  he  mapped  onto  the  plane  tangent  at 
elevation  $0  and  azimuth  A=0.  Hence  the  series  of  points,  specified  by  their 
azimuth  and  elevation  values,  at  which  optical  data  exist  can  be  plotted  on  a  plane 
representing  the  plane  of  tangency. 

To  compare  the  projected  scan  line  with  a  photograph  of  the  cloud  it  is 
necessary  to  know  the  azimuth  and  elevation  of  the  point  of  tangency  and  the  scale  size 
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(i.  e„  the  value  of  R  In  the  above  equations).  From  this  information  cloud  photographs 
should  be  able  to  be  reproduced  which  can  be  directly  compared  with  the  scan  lines. 

It  will  be  noticed  in  the  data  books  that  the  configuration  of  the  scan  grid 
varies  with  reference  elevation.  This  preserves  the  scale  factor,  which  bears  a  direct 
relationship  to  the  magnification  factor,  so  that  a  constant  magnification  can  be  achieved. 
That  if,  a  given  azimuthal  interval  will  take  less  linear  space  at  higher  elevations  than 
at  low. 

The  scan  line  is  presented  at  several  scale  factors,  depending  on  the  angular 
extent  of  the  scan.  The  grid  sizes  used  are:  6°  azimuth  x  6°  elevation;  10°  x  10°,  and 
14°  x  14°.  As  plotted,  the  scale  factors  are: 

6x6  grid;  1”  *  1°  subtense 

10  x  10  grid;  1"  «  1.  67  subtense 

14  x  14  grid;  1”  =  2.33  subtense 
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IQ.  IMAGE  ORTHICON  RECORDS 


Photographs  of  18  frames  of  the  10  record  on  event  Spruce  follow. 

Although  the  quality  of  the  originals  has  been  degraded  through  the  reproduction 
process,  the  striated  structure  of  the  developed  cloud  is  apparent.  The  first  line 
of  the  data  block  contains  the  time  in  hours,  minutes,  and  seconds.  These  records 
can  be  compared  with  the  digital  records  on  Spruce  although  no  attempt  has  been  made 
for  this  summary  report  to  render  the  scale  accurately.  Densitometric  studies  on 
these  photographs  are  very  difficult  because  of  the  non-linear  transfer  function  involved 
in  their  reproduction.  However,  point  values  of  intensity  can  be  deduced  from  intensities 
on  the  video  tape. 

IV.  FORM  OF  THE  REDUCED  DATA 

As  a  sample,  the  reduced  data  from  record  101  of  event  Spruce  is  included 
as  Appendix  2.  Although  in  the  appendix  the  several  frames  of  the  record  are  separate 
sheets,  the  data  books  contain  each  record  as  an  integral  folded  sheet  which  can  be 
removed  for  detailed  work.  This  particular  record  is  shown  as  representative  of  many 
of  the  features  of  the  reduced  data  which  may  be  of  interest. 

First,  the  scan  line  is  plotted  on  a  6°  x  6°  grid.  On  this  record  the  scan 
line  extends  beyond  the  acceptable  limits  for  the  6°  x  6°  grid  thus  it  is  next  plotted  on 
the  10°  x  10°  grid.  The  beginning  of  the  scan  is  denoted  by  a  square.  Each  75th  point 
is  overstruck  with  an  X.  The  data  block  contains  the  event  name,  record  number,  and 
time  at  the  beginning  and  end  of  the  record.  It  also  contains  other  features  which  need 
some  explanation.  The  reference  coordinates  given  are  the  azimuth  and  elevation  of  the 
record  point  at  the  center  of  the  scan  line  rounded  off  to  the  nearest  degree.  The  point 
count  contains  the  total  number  of  points  in  the  record  and  the  number  of  points  plotted. 
The  number  of  points  plotted  is  always  smaller  than  the  total  number  because:  a)  the 
first  and  last  points  are  always  discarded  to  avoid  a  commonly  observed  spurious 
deviation  and  b)  those  points  for  which  the  pressure  setting  of  the  FP  interferometer 
deviate  significantly  from  the  mean  pressure  for  the  record  are  discarded.  The 
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mean  pressure  is  determined  by  averaging  the  last  ten  points  of  the  record  left  after 
the  very  last  point  is  discarded. 

On  succeeding  frames  the  measured  data  is  plotted  against  path  distance 
from  the  initial  point  measured  in  subtended  degrees  along  the  path.  This  has  the 
effect  of  straightening  out  a  linear  scan.  A  specific  point  can  be  located  roughly  by 
estimating  its  position  between  two  X  marks  or  accurately  by  counting  the  points 
from  the  nearest  X  mark. 

In  the  integral  folded  version  of  the  record  the  data  from  specific  channels 
are  presentee1,  together  so  that,  for  example,  the  data  from  the  4934A  channel  of  the 
FP  are  along  the  top  half  of  the  first  part  of  the  record.  This  means  that  in  the  appendix 
these  data  appear  on  the  top  half  of  each  of  the  first  4  pages  of  the  data  presentation. 

Data  are  presented  for  the  FP  4934  channel  the  FP  6142  channel,  the 
photometer  channel,  and  for  the  point  by  point  ratio  of  4934  to  photometer  intensities. 

The  left  hand  axis  of  each  plot  is  the  intensity  relative  to  the  maximum 
intensity  in  that  channel  for  the  specific  record.  This  has  the  effect  of  spreading  the 
data  more  nearly  over  the  entire  useable  area.  The  right  hand  axis  is  the  absolute 
intensity  as  derived  from  calibration.  The  calibration  procedure  is  described  in 
detail  below.  The  ratio  of  FP  4934  to  the  photometer  intensity  has  no  readily  useable 
significance  by  itself  but  variations  in  this  ratio  over  the  scan  are  of  some  interest. 

For  this  reason  the  scale  applied  has  no  quantitative  significance  but  is  effectively 
arbitrary. 

i 

For  each  channel  the  maximum  number  of  counts  is  printed  out.  This  is 
the  number  of  counts  associated  with  the  highest  intensity  point.  Its  main  significance 
is  to  give  an  indication  of  the  statistical  (  fluctuation  expected. 

The  data  plotted  have  been  reduced  by  the  channel  dark  count  determined 
during  the  specific  event.  For  the  second  FP  channel  and  for  the  late  cloud  developments 
in  all  channels,  the  data  count  ratio  is  not  very  much  greater  than  the  dark  count  rate. 
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Hence  the  significance  of  the  data  under  such  conditions  is  highly  questionable. 

V.  APPLICATION  OF  THE  REDUCED  DATA 

Many  uses  of  the  reduced  data  have  been  planned  and  others  suggest 
themselves.  The  primary  use  will  probably  be  a  comparison  of  interferometer  data 
with  intensities  derived  densitometrically  from  photographs  of  the  cloud.  To  accom¬ 
plish  such  a  comparison  it  is  necessary  to  determine  the  points  on  the  photograph  to 
which  the  interferometer  data  apply.  This  can  be  done  by  using  the  overlays  supplied 
in  conjunction  with  photographs  reproduced  at  the  appropriate  scale. 

In  many  cases  a  series  of  spatial  scans  covers  a  nearly  constant  portion 
of  the  cloud  and  such  a  series  can  be  useful  in  performing  frame  to  frame  correlations 
of  structure.  A  preliminary  examination  of  frames  68  to  71  of  event  Spruce,  done 
before  the  reduced  data  were  available  in  the  present  format,  has  yielded  the  qualitative 
conclusion  that  certain  features  do  correlate  from  frame  to  frame  over  a  period  of  a 
few  minutes  while  certain  other  features  may  appear  in  one  or  two  frames  but  shortly 
vanish  from  view.  With  the  data  in  its  present  format  more  detailed  study  of  the 
evolution  of  striae  and  their  short  term  behavior,  at  least  within  the  plane  normal  to 
the  line  of  sight,  is  possible. 

The  cloud  expansion  in  its  early  stages  can  be  studied  and  non-symmetric 
models  for  the  expansion  process  can  be  matched  against  the  observed  intensities  over 
the  various  scans  and  scan  directions  used.  Edge  steepening  can  be  further  studied  by 
using  the  data  plots. 

VI.  CALIBRATION 

The  calibration  in  the  field  was  accomplished  using  both  a  low  brightness 
radioactive  phosphor  (LBS)  which  bad  been  calibrated  against  a  standard  at  ESSA  some 
6  months  before  use  (Reference  ;?  )  and  an  electroluminescent  source  (ELS)  of  con¬ 
siderably  greater  brightness  which  was  calibrated  against  the  LBS  on  each  run. 


7 


Calibration  was  accomplished  at  the  end  of  each  event  (with  the  exception 
of  Plum)  and  consisted  of  taking  dark  counts  on  all  channels,  photometer  response  to 
LBS,  photometer  response  to  ELS,  and  Fabry-Perot  response  to  ELS. 

J'he  LBS  was  not  large  enough  to  fill  the  entire  photometer  aperture,  so 
the  effective  LBS  brightness  for  full  aperture  was  reduced  by  the  ratio  of  the  areas. 

The  LBS  had  a  2. 75  inch  active  diameter,  and  the  photometer  aperture  was  4. 30  inches 
in  diameter,  so  the  ratio  factor  is  2.44. 

The  radiance  of  the  cloud  (above  the  atmosphere)  as  seen  by  the  photometer 
can  then  be  written: 


<NPC  ~  NPP) 

<npl  -  npd> 


2.44  T, 


where  I  is  i.e  LBS  intensity,  and  T.  is  the  atmospheric  transmission  (in  which 

A 

the  effect  of  the  plate-glass  observation  window  are  included),  and  the  subscripts 
on  N  indicate  counting  rates  of  the  photometer  (P)  when  looking  at  the  clouds  (C),  the 
LBS  source  (L)  or  dark  (D).  The  LBS  intensity  incident  on  the  photometer  is  determined 
by  multiplying  the  measured  continuum  radiance  of  the  LBS  per  angstrom,  iL,  by  the 
full-width-at-half-height  of  the  particular  interference  filter,  Wp,  in  angstroms.  The 
intensity  then  becomes 


(NPC  "  NPP) 

<npl  -  npd> 


*L  WF 
2.44  T , 


Intensities  derived  thusly  should  be  directly  comparable  with  (filter)  photographically 
derived  cloud  radiances  if  both  correspond  to  the  same  viewing  aspect  of  the  cloud. 

As  Hake  et  al  have  pointed  out  (Reference  3  ),  calculation  of  the  radiance 
of  the  cloud  from  the  FP  signals  is  a  more  complicated  matter  because  the  continuum 
emission  from  the  ELS  is  further  attenuated  by  the  action  of  the  FP.  For  the  purposes 
of  this  calibration  it  is  assumed  that  the  portion  of  the  continuum  emission  passing 
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the  FP  is  determined  by  the  FP  finesse  F 
the  instrumental  full-width-at-half-heightj 


F  =  free  spectral  range  (FSR)  divided  by  the 
in  the  relation: 


where  ig  and  Ig  are  ELS  intensities.  The  lower  case  ig ,  as  mentioned  above, 
denotes  intensity  in  some  frequency  interval.  In  the  case  above,  the  frequency 
interval  is  one  FP  passband: 

Wpp  =  FSR/F 
or 

wfp<hZ)  «  CWFP 

where  FSR  *  (1/2)  t  is  the  free  spectral  range  of  the  FP,  t  is  the  spacer  length, 
and  c  is  the  speed  of  light.  (The  equation)  is  dimensionally  correct  only  if  it  is 
realized  that  the  ELS  intensity  I_,  is  effectively  distributed  over  one  FSR  of  the 

hi 

FP— i.  e. ,  Ig  ■  intensity/FSR  (Hz). 

F  was  measured  pre  and  post  event  at  the  laser  wavelength,  A. 6328  R, 
and  must  be  computed  for  the  other  wavelengths  using  FP  parameters  (plate  flatness, 
primarily)  established  by  the  A6328  R  meacurement.  If  the  pre  and  post  event  values 
differ  significantly  (>  1  unit),  then  linear  interpolation  should  be  used  to  give  F  for  a 
particular  time.  The  values  determined  for  F  are: 


Wavelength  (R) 

SPRUCE 

F(Wpp/FSR) 

REDWOOD 

F(Wpp/FSR) 

OLIVE 

F(WFp/FSR) 

6328 

17.0 

measured 

17.3 

measured 

19.0 

measured 

6142 

18.4 

computed 

18.9 

computed 

20.3 

computed 

5535 

18.4 

computed 

18.9 

computed 

20.4 

computed 

4934 

16.  6 

computed 

16.0 

computed 

17.7 

computed 
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The  effective  ELS  intensity  for  the  FP  is  determined  by  calibration  against 
the  LBS  using  the  photometer.  The  relationship  is 

,  ‘l  wf  *w-”ra> 

E  -  2.05  x  <Npl  .  NpD) 

The  factor  2. 05  results  from  the  aperture  diameter  of  the  FP  being  100  mm  instead 
of  the  4. 30  inches  appropriate  to  the  photometer. 

The  cloud  radiance  (above  the  atmosphere)  as  determined  by  FP  at  a 
particular  point  in  the  spectral  line  is 

i  =  i  x  ~  NfP).  x  _J _ 

c  E  <nfe-nfd>  tatfp 

where  is  again  the  cloud  radiance  in  a  limited  spectral  width,  the  passband  of  the 
FP.  TFp  is  an  FP  transmission  factor,  that  arises  because  the  FP  was  calibrated 
through  the  wide  field  port  (no  telescope)  but  the  cloud  was  observed  through  the 
telescope.  Tpp  has  been  taken  to  be  0. 74,  from  an  estimate  that  includes  the 
obstruction  by  the  secondary  mirror  and  reflection  losses  from  the  two  additional 
mirrors  in  the  light  path. 

The  dark  counts  have  been  determined  to  be 


EVENT 

PHOTOMETER 

FP  1 

FP  2 

Spruce 

5,  67 

4.91 

2o.  2 

(15) 

Olive 

4.28 

4.62 

21.1 

(15) 

Redwood 

3.09 

3.  58 

17.08 

(10) 

On  FP  2  the  dark  counts  were  so  large  and  the  signals  generally  so  small  that  much 
of  the  data,  after  the  dark  count  is  subtracted,  ends  up  negative  and  thus  does  not 
appear  on  the  data  plots.  To  avoid  this  loss  of  data  points  the  dark  count  actually 
subtracted  is  arbitrarily  chosen  to  be  the  values  in  parentheses  above. 
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The  radiancy  values  have  been  determined  absolutely  according  to  the 
following  relationships. 

Spruce: 

Photometer  I  =  8. 25  x  10-6  N* 

FP  1  I  =  2.47  x  10“15N* 

FP  2  I  =  1. 1  x  10~13  N* 

Redwood: 

Photometer  I  =  8. 45  x  10~3  N* 

FP  1  I  =  2.84x10“15N* 

FP  2  I  «  1. 89  x  10~13  N* 

Olive: 

Photometer  I  *  8.  98  x  10"6  N* 

FP  1  I  -  2.  62  x  10"15  N* 

FP  2  I  «  7. 25  x  10'13  N* 

o 

The  units  are  for  I  (photometer),  erg/cm  sec  sr.  and  for  I  (FP),  erg/cm 
sec  sr.  Hz.  N*  is  the  difference  between  actual  counts  and  dark  counts  each  0. 1  sec 
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APPENDIX  II 


SAMPLE  DATA  SPRUCE  -  #101 


00 


tm 


"I  iinrnin  liMniiMumilf 


RELATIVE  INTENSITY 


INTENSITY  -  ERG/CM* *2-S£C-SR-HZ  INTENSITY  -  ERG/CM 


INTENSITY  RELATIVE  INTENSJ 


INTENSITY  -  ERG/CM* *2-SEC-SR-HZ 


RELATIVE  INTENSITY  RELATIVE  INTENSITY 


INTENSITY  -  £RG/CM**2-SEC-SR-HZ 


V£  INTENSITY  RELATIVE  INTENSITY 


INTENSITY  -  ERG/CM**2-SEC-SR-HZ 


APPENDIX  IH 


PLOTTING  PROGRAM 


Preceding  page  blank 
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PROGRAM  TO  PLOT  SE8BDE  DATA  TAPES 


written  bt  g.lieblJng  -  general  electric  -  aug 


csfcd 

c 

c 


PROGRAM  TO  PLOT  SECEDE  DATA  TAPES 
WRITTEN  BY  G.LIEBLING  -  GENERAL  ELECTRIC  -  AUG  1771 


EXTERNAL  TABLlV 

DIMENSION  T I TlE ( 12 ) #  LABEL(6,4).  NqR  C 100  )  •  IDATA(26), 

PRES {  70  0  1 1  SEC ( 700  J  ,  TmJN(700>, 

A  Z  {  70  0  )  *  EL  ( 700  ) »  FMON(?O0),  FPK700  ).  FP2(?00). 

RAZ  ( 700  I »  REL<700  1  .  RMONC700).  RFPkJOO).  RFP2(700  )',‘ 
RSEC(TOO).  RM I N (  700  ) »  rats;7ooi» 

PARRA Y ( 7flg  ,  4  )  ,  NAME(3) 

DIMENSION  I XTA8( 5g ) ,  IYTaB(30),  DISC  7qq> ,  XTAB(700),  YTAB(700) 
DIMENSION  ASCALE ( 3 )  ,  FMAX<3),  MARTOP<3),  HARBOR).  NTfC<3l, 

1  TVAL(10.  3j.  I L A 8 C & • 3  ) 

equivalence  ( parra y ( i * i ) »  rfpkiJ).  <parrav(i?3) ,  rfpxh), 

1  (PARRAY(1,2),  RMON(l)),  ( PARR  A Y ( 1«  4 ) ,  RATS(1)> 

EQUIVALENCE  ( I  DA  TA  d  >  *  I  DAT  1 J ,  ( I  DA  Ta I  2 ) »  IDAT2), 

1  < I  DAT  A ( 3 ) »  IDAT3),  (  I  DA  TA ( 4  )  ,  IDAT4),  (  J  DATA (5 ) .  I D A  T5  )  » 

2  ( I  DAT  A ( 6 ) ,  I  0  A  T  6 ) f  {  f  D  A  T  A  {  7  )  ,  I  DA  T7 ) ,  {IDATA(8),  I  DATS  ) , 

»TA<11).  IDATui 


DATA 

(( LABEL ( l i J  )  •  IP  1 1 6  )  ,  J  a  1.4) 

1 

/  6  W 

6HFABRV  ,  AHPEROT  , 

6M4934  C* 

2 

6h 

6H  PRO,  6^ TOME  TE  t 

AMR  CHA N t 

3 

6  H 

ahfabry  ,  amp^rot  , 

6M6142  c. 

4 

6HFABRT  ; 

6HPER07  .  6H4934  c. 

AHHANNEL. 

DATA 

( ( I L AR ( I , J ) .  I  *  1,6),  J  «  1,3) 

1 

/6WINTENS. 

6HITV  -  ,  'HERC/CM, 

AW»»2-SE. 

2 

6H!\TENS. 

6  H  I  T  V  -  ,  AMFRG/CM, 

6Wa»?-5E » 

s 

6Wf\TENS. 

6 H  f  Ty  „  ,  6WER0/CM, 

AM»*?~SE  , 

DATA 

RADS/0 ,0174533/ 

n  A  T  A 

IN.  I  OUT,  I  TP  /  5|  6<  l3  / 

scale 

FACTORS  FOR  GRIDS 

SCALFl  =  0,0?0 

DELI 

s  3.0 

DEL  2 

*  5.0 

DEL3 

*  7,0 

MARG 

=  60 

L  I  NC 

s  30 

ldrau 

*  linc  « 

X 

FL  INC 

s  LINC 

ONE  ® 

1,30 

CALL 

C w S 1 7 V  (2, 

2  ) 

CALL 

R I TSTV  (12, 

2d.  7 A8l 1 V  ) 

READ 

UN,  5030) 

title 

VRITFI t  OUT , 6010 I 

Utlf 

READ 

(IN,  5050) 

NAMg 

rf  a  n 

(IN,  5 1 4 0 ) 

[J  A  P  K  r*  j  *  r  R  ’i  ^  f 

read 

(IN,  5 1 6  0  ) 

ASCALE ( i ) ,  A$:alF(3>,  ASCAL! 

6HMANNEL 


Ah 

A* 


6HHANNpL,  6M 

6H/pHoT0,  ftKMgTpR 


6HC-SR  , 


6HZ 

A* 

6HZ 
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C  READ  FILE  NUMBER  OF  DESIRED  FILE 
READ  (IN.*  5100)  NFILE 

C  READ  NUMBER  OF  RECORDS  TO  RE  PROCESSED 
C  AND  INDEX  NUMBER  OF  EACH  RECORD-MAX  I  MUM  100  RECORDS 
READ  (IN,  5120)  NRECS 
READ  (IN,  512q)  (NOR(I),  I  *  i,  NRECS) 

C 

C  SET  Up  LflOp  TO  PROCESS  EACH  RECORD 
READ  (ITP)  IDATA 
DO  750  NR  *  1,  NRECS 
NREC  *  NORtNR) 

GO  TO  205 

C  READ  ID  FRAME 

200  READ  (ITP)  IDATA 

205  IF  (IDAT1  -  97  )  200,  210,  200 

C  CHECK  IF  DESIRED  FILE 

C  CHECK  IF  DESIRED  RECORD 

210  IF  ( I D A T 2  .  NFILE)  200,220,215 

C 

C  FILE  HAS  BEEN  BYpASSED-  ABOrT  JOB 
215  WRITE  ( I  OUT ,  6100)  NFILE,  IDAT2 
REWIND  ITP 
STOP 
C 

220  IF  (IDAT3  -  NrEC)  200.  230,  225 
C  DESIRED  RECORD  HAS  BEEN  PASSED-  TRY  NEXT  RECORD 
225  WRITE  (IOUT.  6050)  NREC 
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PROGRAM  TO  PLOT  SECEDE  DATA  TAPES 
WRtTTEN  B V  G.LIEBL1NG  -  GENERAL  ELECTRIC  -  AUG  197l 


SEC(NF)  «  TEMP2  ♦  (TEMPI  /  100.0) 

TMIN(NF)  a  IDAT1 
300  CONTINUE 

C  IF  THIS  LOOP  COMPLETED  MUST  HaVE  READ  700  POINTS 

c  too  many  -  disregard  Record 
WRITE  C  I  OUT  *  6060)  NREC 
GO  TO  750 

C  INVERT  RECORD  AND  EXTRACT  CONSTANT  PRESSURE  PART 

C  GET  AVERAGE  OF  LAST  TEN  PRESSURE  POINTS 
32(1  Nil  s  NWDS  -  in 
Ml  r  NWDS  o  1 
M2  s  NWDS  -  2 
APRES  =  0.0 
DO  330  I  *  Nil,  N 1 
330  APRES  =  APRES  ♦  PRES  (  l  ) 

APRES  *  APRES  /  10.0 

C  INVERT  RECORD  -  REJECT  FIRST  AND  LAST  POINTS 

NREJ  =  0 

NREV  *  0 

DO  370  IV  a  1.  N2 

1  REV  s  NWDS  -  IV 

IF  (  ABS ( PRES ( I  RE  V )  -  APRES)  ,LE.  4,0)  GO  TO  350 
0  IF  5  CONSECUTIVE  POINTS  REJECTED  DISREGARD  REST  OF  RECORD 

NRFJ  =  NREJ  *  1 
IF  ( NREJ  -  5)  370.  400.  400 
35.)  NREJ  *  0 

NRF V  *  NREV  ♦  t 

RSEC  (NREV)  a  SEC  (IrEV) 

RM’N  (NRkV)  =  TM I N ( IRgV) 

RAZ  (NREV)  *  AZ  ( jRFV) 

REL  (NREV)  s  EL  (IREV) 

RMON  (NREV)  *  FMONUPEV)  -  DARKPH 

RF P i  (NREV)  *  FPl  (IREV)  -  D ARKF 1 

RFP?  (NREV)  *  FP2  (IREV)  -  DARKF 2 

370  CONTINUE 

400  IF  (NREV  .f,T.  5)  GO  TO  405 

write  udut,  6040  )  nrec 

GO  TO  750 

c 

C  NORMALIZE  TABRY-pgpOT  INTENSITIES  AND  MONITOR 

C  FIRST  FIND  MAX  VALUES 

4Q5  FIMA/  a  RFPKl) 

F?M AX  a  REP2 ( i  ) 

FMMAX  =  RMON ( j  ) 

DO  4 1 0  I*  ->,  NREV 

IF  ( P  M  0  N (  I  )  .GT.  FMMaX)  FMMAX  s  RMON(I) 

I F  ( R  F  P 1 ( { )  .GT.  FlMAX)  F1MAX  =  RFPl(l) 

4 1  i)  IF  (  RFP2  (  I  )  .GT.  F?MAX)  F  2  M  A  X  c  RFP2(I) 

FMAX(l)  »  FlMAX 
FMAX(2)  a  FMMAX 
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FMAX<  3 )  *  F2MAX 

ABSOLUTE  INTENSITY  SCALE 
DO  413  I  *  1,3 
ABSMAX  *  FMAX(I)  *  ASCALE(I) 

ITEN  S  I  NT  C  ALOGIO(ABSMAX)  -  1.0) 

UNITS  =  ABSMAX  /  i0.0  »»  ITEN 
SINC  x  2,0 

jp  ( UN i TS  ,LT.  6.0)  S  » NC  «  1.0 
IF  (  ijNITs  .LT.  3.0)  s*NC  8  °*5 
NT  I C< I  )  *  UNITS  /  SINC  *  2.0 
SINC  s  SINC  *  10.0  **  ITEN 
NTICR  a  NTICU  ) 

NT  I  Cl  =  NTICR  -  1 
TVALIi.  I)  a  0.0 
TVALINTICR,  I )  *  ABSMAX 
no  413  j  *  2,  NT  j  Cl 
41?  rVAL(J.I)  *  TVALU-1,  I)  *  SINC 
C 

C  REINVERT  arrays 
NREV2  *  NREV  /  2 
DO  415  I  2  1,  NREV2 
NrW  x  NREV  ♦  1  -  I 
TEMP  -  RAZ(  I  ) 

RAZ(I)  s  RAZ(NFW) 

RAZ(NFW)  s  TEmP 
TEMP  2  REL(  I  ) 

RpL(I)  s  RFL<NpW) 

RpL ( NpW )  x  TpMP 

DO  415  J  *  1,  3 
TEMP  x  PARRA  Y (  I  .  J  ) 

PARRAY(!,'  J)  -  PARRAY(NFW,  J) 

4 1 5  PARRAY<  NFW,  J  >  x  TEMP 

C 

c  now  normalize,  also  get  Ratio  of  monitor  to  fpi 

DO  420  1*1,  NREV 

rATS<  I  )  *  rFPi (I)  /  rMON ( I  ) 

RMON (  I  )  *  RMON (  I  5  /  FMMAX 
HFP1 (  I  )  t  RFPl (  I  )  /  F 1 M  A X 
42n  RFP2<  I  )  *  RFPjni  ,  F2MAX 
C 

c  take  appropriate  average  of  fP2  channel 

NAVG  r  5 

IF  (F2MA t  ,GE.  25.0)  NAVG  s  4 

IF  (F2MAX  , GE ,  5g  .  0  )  N A V G  *  3 

IF  (FpMAX  , GE  .  t  0  0  •  0  7  N A  V  G  =  ? 

IF  (F?MAX  .GE.  200.0)  NAVG  =  1 

jF  (NAVG  ,EQ.  n  GO  To  4P4 
FA  vG  =  NA\/G 

NREVA  =  NREV  -  NAVG  ♦  1 
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DO  422  1  «  1,  NREVA 

SUM  *  0,0 

DO  421  J  *  1,  NAVG 

IJ  *  I  ♦  J  -  1 

421  SUM  *  SUM  *  RFP2( I J) 

422  RFP2< I )  *  SUM  /  FA  VG 
NREVA  b  NREVA  «•  1 

DO  423  I  *  NREVA,  NREV 

423  RFP2<  I  1  *  RFP2U-1) 

424  CONTINUE 
C 

c  ready  to  set  up  az  el  plot  on  »gnomonic»  projection 
c  determine  center  coord 

IPASS  8  1 

NHALF  *  NREV  /  2 

AZCENT  S  A1NT(RAZ(NHALF)  ) 

ELCENT  s  AINT(REL(NHALF>  ) 

IF  ((RAZ(NHALF)  -  AZCENT)  .GT.  0.5)  AZCENT  «  AZCENT  *  ..0 

IF  (<REL(NHALF)  _  ELCENT)  .GT,  0.5)  ELCENT  >  ELCENT  *  1,0 

AZCENR  a  AZCENT  •  RADS 

ELCENR  =  ELCENT  *  rADs 

C 

C  determine  min  and  max  ANGLES  and  deltas 

AZMAX  =  0.0 
A Z M I N  b  0,0 
ELMAX  =  o.o 

Elm  in  *  o.o 

DO  47 o  I  *  1,  NREV 

TEMP  t  RAZ(I)  -  AZCENT 

IF  (TEMP  -  AZMAX  )  430,  45q,  425 

425  AZMAX  s  TEMP 
GO  TO  450 

430  IF  (TEMP  _  AZMIN)  440,  450,  453 

440  AZMjN  s  TEMP 

453  TEMP  =  rEL ( I )  -  ELCENT 

IF  (TEMP  -  ELMAX)  460.  47(1,  455 
455  ELMAX  s  TEMP 
GO  TO  470 

460  ;F  (TEMP  .  ELM  I N )  465.  470  .  470 

465  El^iN  a  TEMP 
4  7  o  CONTINUE 

DELA2  *  AZMAX  -  AZMIN 
DELEL  =  ELMAX  -  ELM  I N 
DELANG  a  DEL AZ 

IF  (DELEL  .GT.  DELaZ)  DELaNG  »  DELEL 
C 

C  FIRST  PLOT  ON  6X6  DEG  SCALP 

SCALF  a  SCALEl 
DEL  s  DELI 

C  INSERT  8LANK  FRAME 
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47*5 

/> 

O 

C 

c 

c 


c 


c 

477 

4  80 
C 

c 

c 


49o 

C 

^  V»  y 

C 


CALL  FRAMEV  (3) 

ELROT  s  ELCENT  -  DEL 
FLTOP  *  ELCENT  ♦  DFL 
AZLFT  a  A7CENT  -  DEL 
AZRGT  a  AZCENT  ♦  DEL 
advance  FILM 
CALL  FRAMEV  <2> 

DSC  ALE  s  DEL  *  SCALE 
ESCALE  s  -  DSCALE 

SCALE  GRID  PAPER-  LEAVE  50  UNIT  MARGINS 
CALL  XSCALV  <  FSCALE.  DSCaLE,  MARG.  MARG) 

CALI.  YSCALV  '  -SCALE.  DSCALE.  MARG.  MARG) 

DRAW  LONfiT.  )E  LINES 
DINC  *  DEL  /  2.0 
AZLINE  -  AZLFT 
DO  480  I*  1.  5 

CALL  CDORn  (XL1.  YL1.  AZCENR,  ELCENR.  AZLINE.  ELROT) 

CALL  COORD  (  XL2,  YL2,  aZcEnr,  ELCENR,  aZL’NE,  ELTOP) 

I  XL  1  =  | XV  <XL1> 

T  XL  2  *  IXV  (XL2) 

JYLl  *  I Y V  (YLl) 

IYL2  *  IYV  <YL2) 

CALL  LINFV  (  I  XL).  *  IYL1.  IXL2,  IYL2) 

DETERMINE  LIMITS'  FOR  USE  IN  INTENSITY  pLOTS 

IF  (I  ,NF.  3)  GO  TO  477 

DFL  Y  =  YL2  -  YLl 

DSCL  =  (2.0  *  DEL)  /  DElY 

LAREL  LINE 

I  XL  1  s  I  XL  1  -  24 

IYL1  =  I  Vi_  1  -  10 

CALL  LABLV  (AZLINE,  IXL1,  IYL1.  6.  1.  4  ) 

AZLINE  s  AZLINE  ♦  DlNC 

DRAW  LATITUDE  LINES 

SINCE  LATITUDE  LINES  ARE  ARCS  MUST  TRACE  OUT  IN  STRAIGHT  SEGMENTS 

AINC  =  DEL  *  2.0  /  Fl  INC 

FI. LINE  s  ELBOT 

no  s-:o  1=1,5 

AZLINE  =  AZLFT 

OO  4$0  J  =  1,  LORAW 

CALL  COORD  (  XL i .  YLj,  AZCENR,  ELCENR,  AZLINE,  ELLINE) 

IXTAR(J)  =  [XV  ( XL  1 ) 
lYTAR(J)  =  IYV  (VM) 

AZLINE  =  AZLINE  *  AlNC 

tarlf  or  coordinates  for  latitude  line  computfd  -  now  draw  linf 

DO  500  J  =  2,  I  Dt?A w 

CALL  LINE V  (  JXT48(J-1),  I v  T  A  B ( J - i ) ,  ixT*0(ji,  i y  T  A  R ( j  )  ) 

LABEL  LAT  LINE 
IXL1  =  IXTAB(l)  -  46 
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CALL  LABLV  (ELLINE.  !XL1,  lYfAi(l).  6,  1,  4  > 

52o  ELLINE  *  ELLINE  ♦  DlNC 

c  write  info  block 

call  dblock  (  name,-  nrec,  aZCEnt,  blgent,  rminjnrev),  psecinrev).* 

1  RMIN(l)  #  RSECtn.  NWDS;  NREV,  ARRES) 

c 

C  GRID  DRAWN  AND  LABELED 

C  CONVERT  POINTS  TO  X  -  Y  COORDINATES  TOR  PLOTTING 
C 

DO  560  I  *  1,  NREV 

560  CALL  COORD  (XTA9CI),  YTAB (I  } »  AXCBNR,  ELCEN*,  RAZi I ) «  RELM)  ) 
c  now  plot 

CALL  APLOTV  <  NREV,  XTAB,  YTAB,  1,  1,  1,  42,  I  ERR ) 

C  PLOT  EVERY  75TH  POINT  WITH  SPECIAL  SYMBOL 

CALL  APLOTV  <  1,  XTAB,  YTAB,  1,  1,  1,  63,  I  ERR  I 
IE  (NREV  , LT •  75)  Go  TO  567 
00  565  I  s  75, nrev, 75 
IX  *  NXV  (  XTAB(  I  )  ) 

IY  *  NYV  (  Y T A R ( I  )  ) 

565  CALL  R  I  TE  2  V  (  IX,*  IY,  ±  0  0  0  »  $0,  *,  i*  »  i  HX ,  I  ERR ) 

C 

C  PLOT  FINISHED  -  REPLOT  IF  NEggSSARY  ON  c°^prESSED  SCaLE 
567  GO  TO  (  570,  590,  590  ),  JPASS 
5 7 f)  IF  (DELANG  ,LE.  (2.0  *  DELI)  >  GO  TO  590 

IP  (DELANG  ,GT.  (2.0  *  DEL2  j  ,  GO  TO  580 

DEL  *  DEL2 
IPASS  s  ? 

GO  TO  475 
58;:  deL  *  DeL3 

IPASS  a  3 
GO  TO  475 
59"  CONTINUE 
r 

c  PLOTTING  OP  G.NDMONIC  PROJECTIONS  COMPLETED 

c 

<*■ 

C  NOW  PLOT  INTENSITIES 

r 

r  GENERATE  TABLE  OP  LINEAR  DISTANCES  PROM  HRST  pOlNT 

DIS(l)  *  0,0 
DO  6PQ  I  =  2, NREV 

DlNC  *  SORT  (  (XTA8<n  -  XTABU-1))  **  2  ♦ 

1  (YTAB(M  -  YTAB  C  I  «i )  )  #«  2  )  #  DSCL 

620  D  I  S (  I  }  s  DISC  I  ~  1 )  ♦  DlNC 
C  NUMBER  OP  FRAMES  /  ARRAY  «,  2  DEG  /  FRAME 
FRM  *  olS(NREV)  /  2,0 

{P  (<PRM  -  AjNT(FRM))  .GT,  0.0)  PRM  s  FRM  ♦  1,0 
NFRAME  s  FRM 
C 

C  MARGINS 


wwwyjp  j  i  twwpww  xmwm* 
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C 


C 


C 


C 

C 


C 


625 

630 

C 

C 


top  half 

M ARTOP ( 1 )  «  1023  -  924 
MARBOTQ)  s  544 
BOTTOM  HALF 

MARTOP ( 2  )  s  1023  -  462 
MARBOT ( 2 )  *  82 
MARLFT  k  66 
MARRGT  »  !023  -  89j 
DO  690  ! T M  x  1,3,2 

?  DEGREES  /  FRaMP 
DO  680  IFRAM  s  i ,  NPRaME 
CALL  FRAMEV  (2) 

FRM  *  IFRam 


DLEFT  s  (FRM  -  lt0)  *  2-0 
DRGHT  a  FRM  •  2.0 

CALL  PR  j  NT V  (33.  33HPATH  DiSTaNCE  -  SUBTENDED  DEGREES,  364,  60  ) 
CALL  PRINTV  (27,  27H I NTENS I  Tv  VS.  PATH  DISTANCE,  388,  4n) 

2  ARRAYS  /  FRAME 
DO  680  I  PL  =  1,2 
DRAW  GRID 

call  setmiv  (  marlft,  marrgt,  marboT(IPl),  martopcipl)) 

CALL  GRjDlV  (  2.  DLEFT,  DRRhT,  q.O.  °^E ,  0.50r  0.20,  0,  0. 

1  1,  1,  4,  3  ) 

LINE  *  MaRBOT(iPL)  ♦  86 

CALL  R!TF2v  <30.  LINE,  1000,  180,  1,  19,  1, 

1  IAHRElATjve  intensity,  ierr  ) 

LINE  x  io44  -  MARTQP(IPL) 

I  GW  x  IFM  *  I  PL  -  t 

CALL  PRINTV  (  36,  LAOELtl,  IGR)  100.  LINE  ) 

If  ( IGR  .ED.  4 )  gD  Tq  630 

CALL  PrIN'TV  <  12,  12RMAX  CO.jNt  x  ,  <00,  LINE  ) 

CALL  LA0LV  (  FMAX(IGR),  500*  LINE,  4,  1,  4  ) 

LINE  *  1008  -  M  A  »  T  0  P  ( I P L ) 


CALL  RlTEjV  (  .000,  LINE,  50.  0«  <  ,  3.,  .,  I L  A  3  ( 
DRaW  ABSOLUTE  INTENSITY  TIC  MaRKS  ‘ 


IGR),  |ERR) 


LSTOP  «  I  X V  (  D R r, H T  ) 
LStA»T  s  lSTOP  -  3 J 


LCHAR  *  LS^QP  *  10 

N  T  J  C  P  =  N  T I C ( IGR) 

DO  625  NT  s  1.  \Tjrp 

POS  x  TVAL ( NT,  iOR)  /  Tval(nT|CR,  IGR) 

IPOS  =  N  V  v ( P  0  S ) 

CALL  LINFV  (  lsurt,  JPOS,  lSt?R,  JPOS  ) 

CALL  LABLV  (,  TVAl(NT,  IQR),  LChA*.  JPCS,  -3,  1,  1  ) 
CONTINUE 


plot  points 

CALL  APLOTV  (NREV,  DIS.  P A RR A v ,1  .  I  OR ) ,  1.  1  .  ♦.  45,  tFRR) 
CALL  A  PL  0  T  V  (  1,  DIS,  RaRRa y ( j, '  jsR),  l . " l , ' i , *  6 3 . ” I  ERR > 

IF  (NREV  .LT.  755  G7  y 0  6«3 


I 


\ 

j 

5 


1 


i 


1 


r>:i 


i 

5 


i 
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DO  670  I  *  75,NREV»75 
IX  s  NXV  <  D I  S  <  I  ) ) 

I Y  s  NYV  (  PARR A  Y ( I , I GR )  ) 

\F  <!#  «EQ.  0  .OR.  jY  ,60.  Q)  GO  TO  670 

CALL  R I T^2V  <  IX,  I Y ,  1000,  $0,  1*  1*  ’1*  1HX,  I  ERR ) 

6  7  ()  CONTINUE 
680  CONTINUE 
690  CONTINUE 
700  CONTINUE 

C  END  LOOP  FOR  ONE  SCAN  FRAME 

WRITE  IIOUT,  6070)  NREC 
750  CONTINUE 


REWIND 

I  TP 

r» 

STOP 

C 

5  05  0 

format 

(12A6) 

5100 

FORMAT 

(  I2> 

51^0 

format 

(1015) 

5140 

FORMAT 

(  3  F 1 0  •  3  ) 

5 1  *  n 

FORMAT 

( 3E l 0 . 3 ) 

<>oi  p 

FORMAT 

( 1  HI »  3oX.  12A6  ,  ) 

6  04  0 

format 

(,  35X  7MREC0RD  , 

13 

.3  0  H 

SKIPPED  -  INSUFFICIENT 

POINTS 

6  0  5  C 

tormat 

< /  35X  7HRECORD  , 

13 

i  ?8H 

SKIPPED  -  NOT 

FOUND  ON 

TAPE  ) 

format 

{,  35X  7HREC0RD  , 

13 

,  26  M 

SKIPPED  -  TOO 

MANY  POINTS  ) 

60  ;c 

format 

(/  35X  7HREC0RD  , 

13 

,  9H 

PLOTTED  ) 

610“' 

r  0  R  m  a  T 

(  ///  6H  PILE  l 2 

1 7 H  BYPASSED  -  FILE 

.  1 2  6H 

FOUND- 

t 

/  12W  JOB  ARO^TFD 

# 

) 

f  1 

r  N  !* 
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CCORD 

C 


SUBROUTINE  to  COMPUTE  GNOMON  I C  PROJECTION  COORDINATES 
FROM  INPUT  AZ  AND  EL 

SUBROUTJNC  COORD  (  X,  V,  AZCENR,  ELCENR,  AZ,  EL) 

DATA  RADS/o, 01745331/ 

AZR  «  AZ  #  RADS 


ELR  *  EL  *  RA[)S 
Sl^L  s  SIN(AZR  -  AZCENR) 
COSL  *  COS ( AzR  -  AZCENR) 
SINP  x  SINiELR) 

COSP  x  COS ( ELR  ) 

SINPO  x  SIN(ELCENR) 

COSPO  e  COS(ELCENR) 


DENOM  x  SINPO  •  SINP  ♦  COSPO*  COSP  #  COSL 
X  s  <  SINL  *  COSP)/  DENOM 

Y  s  (  COSPO  •  SINP  -  SINPO  *  qOSL  «  COSP)  /  DENOM 

RETURm 

END 
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01-05-72  13.751  SUBROUTINE  TO  PRINT  HOUSEKEEPING  DATA  BLOCK 

ON  GNOMON  I C  PROJECTION  PLOTS 

CDBLK  subroutine  TO  print  housekeeping  data  BLOCK 

ON  GNOMONIC  PROJECTION  PLOTS 
USES  VECTOR  CHARACTERS  for  larger  letters 
SUBROUTINE  DRLOCK  (NAME,  NREC,  AZceNT,  ELCENT,  RMSRT,  RSSRT 
1  RMSTP,  RSSTP,  NWDS,  NREV,  APRES} 

EXTERNAL  TABLlV 
DIMENSION  NAME ( 3 ) 

ILF  =  790 
IH  x  98fl 
I  DEL*  20 
M4R  =  1023 
CalL  CHSlZV  <2,  2) 

CALL  RjTSTV  (  12.  20.  TABlIV) 


CALL 

R I TE2 V  ( 

ILF, 

IH,  MAR, 

9o. 

1, 

16,  1. 

NAME,  IERR) 

IH  = 

IH  -  IDEL 

CALL 

RITE2V  ( 

ILF, 

IH,  MAR, 

90, 

1. 

13,  -1, 

13HREC0RD  NO.  - 

TEMP 

=  NREC 

CALL 

BNBCDV  ( 

TEMP 

,  BCD,  NDS, 

CALL 

RITE2V  ( 

946, 

IH,  MAR, 

*0, 

1, 

NDS ,  1, 

BCD,  IERR) 

IH  = 

IH  -  1  DEL 

CALL 

RITE2V  ( 

ILF, 

IH,  MAR, 

90, 

1  • 

11*  "1* 

iiHREF  COORD  -  , 

IH  = 

IH  .  I del 

£  ALL 

R  I  T  E  2  V  ( 

ILF, 

IH,  MAR, 

90, 

1, 

5 ,  -1 , 

5«  AZ  ,  I  ERR ) 

CALL 

B.jUCUV  ( 

AZCE 

nT,  BCD, 

NnS) 

CALL 

RITE2V  ( 

850, 

IH,  MAR, 

90. 

1, 

NDS,  1, 

BCD,  IERR) 

IH  = 

IH  -  I  del 

CAlL 

RITF2V  ( 

ILF, 

IH,  MAR, 

90. 

1 , 

5,  -1, 

5H  EL  ,  IERR) 

CALL 

BNBCDV  ( 

ELCENT,  BCD, 

\ins) 

call 

RITE2V  ( 

850  , 

IH,  MAR, 

V 

1' 

NDS,  1, 

BCD,  IERR) 

IH  5 

IH  =  T  DEL 

CALL 

R  I  T  lr  2  V  ( 

ILF, 

IH,  MAR. 

90, 

1  • 

4  1  ~  , 

4  H  T ! M  E ,  IFRR) 

IH  = 

IH  „  IDEL 

C  A  L  L 

RITF2V  ( 

ILF, 

IH,  vAR, 

90, 

1 , 

19,  -1, 

18H  START  m 

IERR) 


CALL  BNBCDv  <  RMSRT,  BCD,  NOS) 

1ST  *  886 

IF  (NDS  • EO •  t)  1ST  x  696 

CALL  R!TE2V  (  1ST.  IH,  MAR,  <*0,  1*  NDS,  1,  BCD,  IERR) 
TEMP  s  aINT  (RSSRT) 

CALL  BNBCDV  (  TEMP ,  RCD,  NDS} 

1ST  =  934 

ir  (NDS  , EO .  1>  1ST  i  946 


CALL 

R  I  T  f  2  V  ( 

1ST, 

IH,  MAR ,  90, 

1, 

NDS,  1, 

BCD,  IERR) 

TEMP 

x  RSSRT  - 

TEMP 

CALL 

BNBCDV  ( 

TEMP  , 

BCD,  \DS ) 

call 

RITE2V  ( 

970  , 

IH,  MAR,  90, 

1  , 

2,  1, 

BCD,  IERR) 

,K  x 

I H  0  j  dfl 

CALL 

p  I T  f  2  y  < 

!  Lr  1 

!  H ,  M  A  R ,  9  0, 

4 

J.  9 

18,  - 1  » 

ion  STOP 

1 

IERR) 

CALL 

BNBCDV  ( 

RMSTP 

,  8CD ,  NDS) 

i  01-05-72  13.751  SUBROUTINE  TO  PRINT  HOUSEKEEPING  DATA  BLOCK 

ON  GNOMON J C  PROJECTION  PLOTS 

1ST  *  886 

IP  ( NQS  ,EQ.  1ST  *  898 

TEH'  J'mnT  (RsItPI  IK'  MAR’  ,0'  NDS’  11  8CD’  1ERR) 

CALL  BN8CDV  <  TEMP,  RCD,  NDS) 

1ST  *  934 

IP  (NDS  ,EQ.  *)  1ST  *  946 

CALL  RIT82V  (  1ST,  IH,  MAR,  *0,  1#  NDS,  1#  BCD,  IERR) 

TEMP  s  RSSTP  -  TEMP 

CALL  9N8CDV  (  TEMP,  BCD,  NDS) 

CALL  RIT62V  (  970,  IH,  MAR,  90,  1,  2,  1,  BCD,  IErrj 

!h  *  [H  «  {DEL 

CALL  ft!TE2v  (  ILF,  IH,  MAR,  90,  1,  13,  -1,  13HP0INT  CO,jNT  -.IERR 

im  c  j  h  «,  I  DEL 

TPMP  - 1 NWDS  (  ILF'  IH’  MAR'  1#  ~1'  7H  T0TAL'  IERR> 

CALL  BNBCDV  (  TEMP,  BCD,  NDS) 

IN1*  ?HTE2 1 pIl 92? *  IW'  MAR'  9°*  1#  NDS>  8CD'  IERR) 

TEMP  *!NR£V  (  lLr'  MAR'  9°‘  1#  *'  m1’  9H  JOTTED.  IERR) 

CALL  BNBCDV  (  TEMP,  BCD,  NDS) 

Vnc.  **22'  lM*  MAR'  9°*  1(  NDS*  1*  BCD ,  IERR) 

In  *  I H  ®  I  DEL 

CALL  RITE2V  (  ILF,  IH,  MAR,  9q,  j,  ,  -l#  ^SHAVERAGE  PRESSURE 

IH  .  IH  .  tDEL  ,6BR1 

CAll  BNBCDV  (  APRES,  BCD.  NDS) 

returnITE2V  f  826'  1W'  MAR'  9o’  1#  N°s'  BCD’  IERR) 

END 
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